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Abstract 

A simple topographic slope correction has been developed using a linear- 
ized thermal model and assuming slopes less than about twenty degrees. The 
correction can be used to analyze individual thermal images or composite 
products such as temperature difference or thermal inertia. Simple curves are 
provided for latitudes of 30 and 50 degrees. The form is easily adapted for 
analysis of HCMM (Heat Capacity Mapping Mission) images using the DMA (Defense 
Mapping Agency) digital terrain data. 

Introduction 

Topographic slope effects commonly have been observed on thermal Infrared 
images and are used to deduce geologic structure and morphology (Rowan and 
others, 1970; Wolfe, 1971; Sabins, 1969; Offield, 1975; Cannon, 1973; 

Schneider and others, 1979). Some early thermal model studies considered the 
general nature of slope effects (Watson, 1970; 1971) and it ms suggested 

(Watsen and others, 1971) that a simple correction to thermal models could be 
applied by representing the north-south slope component as an equivalent 
change in latitude and the east-west component as an equivalent shift in local 
time. 

The initial thermal-inertia study (Pohn and others, 1974), using low- 
resolution meteorological satellite data, assumed that at that small scale 
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■lopes, «nd hence topographic effects, w>re negligible. The first topographi- 
cally corrected thermal Inertia Image was constructed from large-scale air- 
craft data using a lookup table for temperature effects as a function of slope 
and azimuth derived from a thermal model (Gillespie and Kahle, 1977). 

The purpose of this paper Is to present an analytical expression for the 
topographic slope correction &:o thermal Images. This expression will be 
primarily for the Interpretation of data from the HCMM satellite but also 
applicable to aircraft data under appropriate slope conditions. 

Method 

A topographic correction term Is developed using the linearized thermal 
model proposed by Watson (1975) and expanded In later studies (Hiller and 
Watson, 1977; Price, 1977; Pratt and Ellyett, 1979; Watson, 1979). The effect 
of slope Is modeled by assuming that It Is due solely to the modulation of the 
incoming direct solar radiation. Second order effects of diffuse and sky 
radiation, sensible and latent heat, geothermal effects and reradlatlon from 
adjacent surfaces have been Ignored. These assumptions restrict the model to 
small slopes (probably less then twenty degrees) and, although generally valid 
for the 500-m HCMM data, do not apply to large-scale aircraft data acquired in 
areas of rugged relief. 

Let v(x,t) ■ Vjj(x,t) + Vj(x,t) where v is the temperature distribution at 
a depth X below the Earth’s surface and at time t measured from local solar 
noon, and v(x,t) also satisfies the heat conduction equation In one dimen- 
sion. Vq Is the temperature for a flat surface and v^ for the additional 
effect due to an Inclined surface of slope d and direction of slope of slope ^ 
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(measured counterclockwise from north)* The surface boundary condition satis 
fied by Vj is: 


3 V ^ 

■*^Tx " ^dir *d “ *o^ Vj «t X ■ 0 


( 1 ) 


and the one dimensional conduction equation is 


2 

3 V 


3x' 


I 

K 



( 2 ) 


where K and k are the thermal conductivity and diffusivity, A and c are the 
surface albedo and emissivity, is the mean surface temperature for the flat 
surface, a is the Stefan-Boltzmann constant, is the amplitude of the 

incident solar flux (equal to S^t where Sj, is the solar constant and x the 
atmospheric transmission factor), and zj are the zenith angles of the flat 
and inclined surfaces, respectively. Using a Fourier series solution (Carslaw 
and Jaeger, p. 74, 1959) the surface temperature of the inclined surface is: 

_ 3 j B h cos(u)St-6 -0 ) 

Vd(0,t) - (1-A) (4 g07^ )’ I — g (3) 

S"0 2 2 

/ (h+W ) + W 

s s 

where the Fourier coefficients Bg and 6g are derived in Appendix One and the 
radiation constant h, the wave number Wg and the phase shift 6g are defined as 
follows: 

h ■ Aeov /K 
o 
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(Watton, 1979)* 


- ✓ w»/2< 

6^m tan’^(W^/(h+W^)) 

Now because Wg/h>l and By ■ 0 for k ■ 3, 5, ••• then 6, * n/4 and 

» 

cos (wt-0i-n/4) 

v.(o,t) •» Q.. (1-A) (4cov ) { B ; 

^ 6 o ^ ^2 (Wi/h) 

• B-^ cos (2r(i>t-e,^-n/4) 

+ I JL— ^ } 

r-1 / 2 (W 2 j/h) 

Now W 2 y/h ■ / 2r (W^/h) 

_ 3 

where W^/h ■ P / w/2 / (4cov^ ) 
and 6y ■ tan"^ (k tan^/sinX) 

Therefore, v^(0,t) » (l-A)d (4eov^ ) { - 

2 2 — TT 

+ (h/W^ ) [ £9^ cos (ut tan” (tan4>/sinX) - tt/4) 

2/2 

^ (-1)^^ / (2r slndi) ^ (slnX cosj.) 

r-1 s l(2r) - 1) / r 

• cos (2ru)t + tan“\2r tan^i/sinX) - s/4) ] } 
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Introduce 


vjDAY , vj(0,0) 

- Vj(0,T/2) 

«nd Avj - vj(0,0) - vj(0,T/2) 

3 1 

Then » Qj^^(l-A) d (Acov^ )” - cos^. ainX/ir 


•¥ (h/Wi ) [ + I (-1)*^^ (2r iln» + sinX cosit) ] 


r«l 


ir / 2r [(2r) - Ij 


^^NITE ^ (1-A) d (4cov^ ) | - cos^ sinX/w + 


- 3.-1 

t 

0 


rfl 


<h/W ) [ “ cosji slnX ■» sln(t> ^ (-1) (2r slnijt *f slnX cos^) j | 


r-1 


It rir K2r) - 1] 


and iv_ . Q^i,(l-A) d (Aeo^/)'' [ (h/W,) ( } ] 


The temperature-difference effect Is proportional to slope and the azi- 
muth variation, which is the term in brackets ([ ]), and is plotted for the 
day, night, and tempr>'ature difference in Figure 1 using an assumed value of 
(h/Wj). The amplitudes of the three expressions are roughly similar, and this 
shows tViat the temperature-difference expression does not provide any reduc- 
tion of the topographic slope effect as an Improvement over the individual day 
and night Images. There are orientation effects, however, that can cause the 
topographic effect to be more pronounced in our image. An HCKM test site in 
Cabeza Prleta, Arlz. (latitude 30°) has a general NW-SE topographic grain — 
corresponding to ^ • 135° and 315o, xt these azimuths, the topographic effect 
should be least in the temperature-difference Image (and hence thermal Inertia 
which is derived from it) and greatest in the daytime thermal image. 



Conclusion 


The general expression given for the topographic correction^ equation 5, 
can be easily aodlfled for other data acquisition tlaes and latitudes. A.n 
expression for can be obtained for clear sky conditions (Huoner-Mlljer , 

1961) and the value of Wj/h conputed as a Becond''order correction for each 
Image pixel from the teaperature difference and sum. The method Is United to 
areas of low topographic slope— probably 20 degrees or less— and Is most 
applicable to HCMM data that have a 500-n resolution and can be aerged with 
the DMA digital terrain data. The technique Is applicable to aircraft data, 
exclusive of areas with rugged relief, but requires digital terrain data at 
appropriately higher resolution and Is thus not generally available. 
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Figure l.—The azimuth variation of the topographic correction factor is 
plotted versus azimuth for the temperature day and night and tem- 
piirature difference at two latitudes (30° and 50°), assuming a 
value of h/Wj ■ 138 (using c ■ 1, “ 280 K, P “ 1500). 
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Appendix One: Fourier coefficients for the topographic slopt'. effect* 


Let’ f(t) • 


T B. cos (wkt - 6|,) 

k-0 


( 1 ) 


where f(t) 

< cos > 


< cos zj > - < eoB &Q> 

CnS Z^ tpjp<t<tgg 

0 otherwise 


< cos *0 ^ 

■ 0 otherwise 

cos Zj A cosut 't' B sinut + C 
cos z ■ cosX cos6 cosoit <f sinX sin6 

A ■ cosd cosX cos£ - sind cos5 cos^ sinX 
B ■ sind sln4» cos6 

C ■ cosd slnX sin6 + sind cos^ sin5 cos.\ 

(Vatson, 1975;1979) 

ujtj • k/2 + sin”^tanX tan6) 

u>t-..’“ x/2 + MIN(sln“'^(tanX tan6 ) , sin”Uc/R) + tan^^(d/A)) 

99 

ojtgr ■ + MAX(-sln"ktanX tand), -Sin’’kc/R) + tan”HB/A)) 

— Z 7 

and R • ^ A > B 

Because the solar declination is a small angle, Id I v 23.5° and thv^. topograph- 
ic slopes are assumed tc be small |d{ < 20° %ie shall introduce the small angle 
approximations for both 5 and d: 

sin X ■> X 
cos X S 1 
x^ - 0. 
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ThuK 

A ■ cc<X - d coi^ ilnX 
B « d sln^ 

C * ts iinX 
wtj ■ w/2 +6 t«nX 
ojtgg * wti + MIN (0, d iln^/coEX) 
WAX «ln^/cosX) 


From equation 1 and using the standard Fourier coefficient method 


- 2 K I 2 

B. cose. / cos u.kt dt - 2 / cos kx dr. 


T/2 


-T/2 


tii 


0 


t ‘•1 

SB 

- / COS a . cos wkt dt - J cos z cos u»kt dt 

'.r 

“'.s “'1 

■ A / cos X cos kx dx - — ^ / cos x cos kx dx 

u,t -wtl 

" sr 

Wt Wtgg 

+ 1 r sin X cos kx dx + ^ / cos kx dx - sinX-S / cos kx dx 

W i “ ijt -li)tl 


x/2 

, _ 2d eositi sinX r ^ ^,^jg 
u) X 


In a similar fashion 


2B. sine^ IT 2 
— = / iin kx dx 


2d sin i|i 
“ ” 0 


sin X sin kx dx 
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Now / cot X cot kx dx 
0 


■ x/A 


- (-1)'+' / U2r)2-H 


and jf tin x tin kx dx 0 
0 

- x/A 


2r(-l)f+i / I(2r)2-H 


Therefore 


Bq ■ -d cot4» tinX/x 


0n - 0 


B^ ■ (d/2) / tin ^ + cos ^ sin \ 


■ tan"^ (-tan ^/sinX) 


y_i%r+l 


B . 2d (-1) 
®2r — T2 


X l(2r) -U 


/ (2r sln(|i) -¥ (sinX cosi}>) 


®2r ” ^■■2r tan^/slnX) 

®2r+l ® 


IT " Ip 2 p 3 p t • t 


fc* ^ 


Reftrcncet 


Cannon^ P. L. , 1973, The application of radar and Infrarad laagary to quantl*- 
tativc gcomorphlc Investigations: Proc. of Renote Sanalng of Earth Res. 
Second Annual Gonf., Tallahoaui, Tann. , p. 503-526, The Unlv. of Tenn. 
Space Inst., Tullahona. 

Carslaw, H. S* , and Jaeger, J. E. , 1959, Conduction of heat In solids: Oxford 

at the Clarendon Press, Second Edition, 510 p. 

Gillespie, A. R. , and Rahle, A. B., 1977, Construction and Interpretation of a 
digital theraal Inertia laage: Photograounetrlc Eng. and Renote Sensing, 

V. 43, p. 983-1000. 

Hunner-Mlller, Susanne, 981, Slnple expressions for the diurnal variation of 
four ataospherlc paraneters for use In surface tenperature models : 
Submitted to Journal of Geophysical Research. 

Miller, S. H. , and Wataon, Kenneth, 1977, Evaluation of algorithms for geolog- 
ical thermal-inertia napping: Proceedings 11th International Symposium 

Remote Sensing of Environment, v. 2, p. 1147-1160. 

Offleld, T. W. , 1975, Thermal-infrared Images as a basis for structure map- 
ping, front range and adjacent plains in Colorado: Geological Society of 

America Bulletin, v. 86, p. 495-502. 

Pohn, H. A., Offleld, T. W., and Watson, Kenneth, 1974, Thermal-Inertia map- 
ping from satellite discrimination of geologic units In Oman: U.S. 

Geological Survey Journal of Research 2, no. 2, p. 147-158. 

Pratt, D. A., and Ellyett, C. D. , 1979, The thermal Inertia approach to map- 
ping of soil moisture and geology: Remote Sensing Env. 8, p. 151-168. 


11 


Price, J. C. , 1977, Ihcmal InertU Mpplng: A new view of the Earth; 

Journal Geophyelcal Ret* 82, p* 2582-2590. 

Rowan, L* C. , Of field, T. U* , Vlatton, Kenneth, Cannon, P. J*, and Uttaon, R. 
D. , 1970, Therauil infrared inveatigationa, Arbuckle Mountains, 

Oklahoma: Geological Society of America Bulletin, v, 81, p. 35A9-3562. 

Sabina, F., 1969, Thermal Infrared imagery and ita application to atructural 
mapping in Southern California; Geological Society of America Bulletin, 
V. 80, p. 397-A04. 

Schneider, S. R. , McGinnia, D. F. , Jr., and Pritchard, J. A., 1979, Use of 
aatellite Infrared data for geomophology atudiea: Remote Sensing Envi- 

ronment, V. 6 , p. 313-330. 

Watson, Kenneth, 1970, A. Introduction and auumary, B. Data analysis tech- 
niques, Part 1, Remote sensor application studies progress report, 1968- 
69; Clearing house for Fed. Scl. and Tech. Info (CFSTI), Doc. PB-197- 

098. 

Watson, Kenneth, 1971, A computer program of thermal modeling for interpreta- 
tion of Infrared images: NTIS PB-203-578. 

Watson, Kenneth, 1975, Geologic applications of thermal infrared images; 
Proceedings IEEE 63, no. 1, p. 128-137. 

Watson, Kenneth, 1979, Thermal phenomena and energy exchange in the environ- 
ment, in Mathematical and Physical Principles of Remote Sensing: Centre 

National D* Etudes Spatlalea, Toulouse, France, p. 10A-17A. 

Watson, Kenneth, Rowan, L* C. , and Of field, T. W. , 19)1, Applications of 

thermal modeling In the geologic Interpretation of IR images: Proceed- 

ings Seventh International Symposium Remote Sensing of Environment, Ann 
Arbor, Michigan, v. 3, 2017-20A1. 


Uolfe, Et V. y 1971, Themal IR for geology: Photograflunetric Eng., v. 37, p. 

43-52. 


13 



